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Abstract 

The  flow  field  generated  by  a  transom-stern  hullform  is  a  com¬ 
plex,  broad-banded,  three-dimensional  phenomenon  marked  by  a  large 
breaking  wave.  This  unsteady  multiphase  turbulent  flow  feature  is  dif¬ 
ficult  to  study  experimentally  and  simulate  numerically.  The  results  of 
a  set  of  numerical  simulations,  which  use  the  Numerical  Flow  Analysis 
(NFA)  code,  of  the  flow  around  the  Model  5673  transom  stern  at  speeds 
covering  both  wet-  and  dry-transom  operating  conditions  are  shown  in 
the  accompanying  fluid  dynamics  video.  The  numerical  predictions  for 
wet-transom  and  dry-transom  conditions  are  presented  to  demonstrate 
the  current  state  of  the  art  in  the  simulation  of  ship  generated  breaking 
waves.  The  interested  reader  is  referred  to  Drazen  et  al.  (2010)  for  a 
detailed  and  comprehensive  comparison  with  experiments  conducted 
at  the  Naval  Surface  Warfare  Center  Carderock  Division  (NSWCCD). 


1  Computational  Method 

The  Numerical  Flow  Analysis  (NFA)  code  provides  turnkey  capabilities  to 
model  breaking  waves  around  a  ship,  including  both  plunging  and  spilling 
breaking  waves,  the  formation  of  spray,  and  the  entrainment  of  air.  A  de¬ 
scription  of  NFA  and  its  current  capabilities  can  be  found  in  Dommermuth 
et  al.  (2007);  O’Shea  et  al.  (2008);  and  Brucker  et  al.  (2010).  NFA  solves 
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the  Navier-Stokes  equations  utilizing  a  Cartesian-grid  formulation.  The  flow 
in  the  air  and  water  is  modeled,  and  as  a  result,  NFA  can  directly  model 
air  entrainment  and  the  generation  of  droplets.  The  interface  capturing  of 
the  free  surface  uses  a  second-order  accurate,  volume-of-fluid  technique.  A 
cut-cell  method  is  used  to  enforce  no-flux  boundary  conditions  on  the  hull. 
A  boundary-layer  model  has  been  developed  (Rottman  et  al.  2010),  but  it 
is  not  used  in  these  numerical  simulations,  and  as  a  result,  the  tangential 
velocities  are  free  to  slip  over  the  hull.  NFA  uses  an  implicit  sub-grid  scale 
(SGS)  model  that  is  built  into  the  treatment  of  the  convective  terms  in  the 
momentum  equations  (Brucker  et  al.  2010).  A  surface  representation  of  the 
ship  hull  is  all  that  is  required  as  input  in  terms  of  hull  geometry.  The  nu¬ 
merical  scheme  is  implemented  in  a  distributed  memory  parallel  computing 
environment  using  Fortran  90  and  the  Message  Passing  Interface  II  (MPI2). 
Relative  to  methods  that  use  a  body-fitted  grid,  the  potential  advantages  of 
NFA’s  approach  are  significantly  simplified  griding  requirements  and  greatly 
improved  numerical  stability  due  to  the  highly  structured  grid. 

2  Domain,  Grids,  Boundary  and  Simulation  Con¬ 
ditions 

The  Model  5673  transom  is  shown  in  figure  ,  and  the  relevant  parameters, 
including  the  length  of  the  model,  the  depth  of  the  transom,  and  the  weight, 
are  given  in  Table  1. 

Table  provides  details  of  the  transom-stern  simulations,  including  the 
speed  of  the  model,  the  Froude  and  Reynolds  numbers,  as  well  as  the  sinkage 
and  trim.  For  3.60  m/s  (7  knots),  the  transom  is  partially  wet,  and  for 
4.12  m/s  (8  knots),  the  transom  is  dry.  All  length  and  velocity  scales  are 
respectively  normalized  by  the  model’s  length  (Lo)  and  speed  (Uo).  The 
number  of  grid  cells  along  the  x,  y,  and  z-axes  are  respectively  denoted  by 
Nx;  Ny ;  and  Nz.  The  number  of  sub  domains  and  processors  along  the  x,  y, 
and  z-axes  are  respectively  denoted  by  nf,  rij  ;  and  nk ■  For  the  simulations 
shown  in  the  fluid  dynamics  video,  discussed  herein,  Nx,  Ny,  Nz  =  2688, 
1024,  384  and  n*,  rij,  nk  =  21,  8,  6. 

The  width,  depth,  and  height  of  the  computational  domains  are  respec¬ 
tively  6.0,  1.6983,  0.66667,  ship  lengths  (Lo).  These  dimensions  match  the 
cross  section  of  the  NSWCCD  towing  tank  (Drazen  et  al.  2010).  The  com¬ 
putational  domain  extends  4  ship  lengths  behind  the  transom  and  1  ship 
length  ahead  of  the  bow.  The  fore  perpendicular  and  transom  are  respec¬ 
tively  located  at  x  =  0  and  x  =  —  1.  The  z-axis  is  positive  up  with  the  mean 
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Figure  1:  Model  5673.  Rendered  at  scale.  Grid  line  separation  is  1.524  m 
(5  ft). 
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Table  1:  Table  1  Model  5673  Details. 


Length  Overall,  Lo 

9.144  m  (30  ft) 

Waterline  Length 

9.144  m  (30  ft) 

Extreme  Beam 

1.524  m  (5  ft) 

Bow  Draft 

Fixed 

Stern  Draft,  Tmax 

Variable 

Construction 

Fiberglass 

Displacement 

771.1  kg  (1700  lb) 

waterline  located  at  z  =  0.  A  plane  of  symmetry  is  not  used  on  the  centerline 
of  the  hull  because  small-scale  turbulent  structures  are  adversely  affected. 
Grid  stretching  is  employed  in  all  directions.  Details  of  the  grid-stretching 
algorithm  are  provided  in  Dommermuth  et  al.  (2007).  The  smallest  grid 
spacing  is  0.0005  near  the  hull  and  mean  waterline,  and  the  largest  grid 
spacing  is  0.01  in  the  far  field.  The  numerical  simulations  are  slowly  ramped 
up  to  full  speed.  The  period  of  adjustment  is  To  =  0.5  (Dommermuth  et 
al.  2007).  Mass  conservation  is  ensured  using  the  regridding  algorithm  that 
is  implemented  by  Dommermuth  et  al.  (2007).  Density- weighted  velocity 
smoothing  is  used  every  400  time  steps  using  a  3-  point  filter  (1/4,  1/2,  1/4) 
(Brucker  et  al.  2010).  The  nondimensional  time  step  is  t  =  0.00025. 

The  simulations  are  run  for  30,000  time  steps,  corresponding  to  7.5  ship 
lengths,  on  the  SGI  Altix  ICE  at  the  U.S.  Army  Engineer  Research  and 
Development  Center  (ERDC).  The  data  sets  are  so  large  that  only  time 
steps  20,0000  through  30,000  are  saved  every  40  time  steps  for  the  purposes 
of  post  processing.  The  1.06  billion  cell  simulation  takes  about  90  hours  of 
wall-clock  time  to  run  30,000  time  steps  using  1008  processors.  The  wall- 
clock  time  can  be  cut  in  half  by  doubling  the  number  of  processors  because 
NFA  scales  linearly.  Alternatively,  increasing  the  number  of  processors  to 
10,000-20,000  will  enable  numerical  simulations  of  breaking  ship  waves  and 
Tsunamis  with  25-50  billion  grid  cells  within  the  next  year. 

3  Flow  description 

The  fluid  dynamics  video  compares  perspective  views  of  laboratory  and  NFA 
results  for  3.60  m/s  (7  knots)  and  4.12  m/s  (8  knots).  The  0.5  isosurface 
of  the  volume  fractions  are  shown  for  the  NFA  predictions.  The  transom  is 
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Table  2:  Calculated  trim  angle,  Froude  numbers  based  on  ship  length,  Fr l 
Reynolds  number  based  on  ship  length,  ReL,  and  transom-condition. 


Speed 

FrL 

ReL 

Trim 

Transom 

(kts) 

(deg) 

(m) 

Condition 

7 

0.38 

3.29  x  107 

0.48 

Wet 

8 

0.43 

3.77  x  107 

0.67 

Dry 

partially  wet  for  3.60  m/s  (7  knots)  and  fully  dry  for  4.12  m/s  (8  knots). 
A  glassy  region  is  evident  behind  the  transom  at  the  4.12  m/s  (8  knots) 
speed.  Significant  air  entrainment  occurs  for  the  3.60  m/s  (7  knot)  case  at 
the  transom,  in  the  rooster-tail  region,  and  along  the  edges  of  the  break¬ 
ing  stern  wave.  For  the  4.12  m/s  (8knot)  case,  air  entrainment  first  occurs 
on  the  forward  face  of  the  rooster  tail  and  along  the  edges  of  the  break¬ 
ing  stern  wave.  For  both  the  3.60  m/s  (7  knots)  and  4.12  m/s  (8  knots) 
speeds,  the  measured  mean  profile  of  the  free-surface  elevation  agrees  well 
with  instantaneous  predictions.  The  fluid  dynamics  video  is  available  at 
www.  saic.  com  /  maritime  /  nfa. 

4  Video 

The  full  size  video,  mpeg4  encoded,  is  approximately  140M6  (download). 
The  web  size  video,  mpeg4  encoded  is  approximately  9 Mb  (download). 


The  videos  are  also  available  at  www.saic.com/maritime/nfa. 

SCENES: 

Fr  =  0.38,  partially  wetted  transom: 


1.  -  Stern-quartering  view. 

2.  -  Stern  view. 
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3.  -  Starboard  view. 


4.  -  Stern-quartering  view  compared  to  snapshot  of  experiments. 

5.  -  Stern  view  looking  down  compared  to  snapshot  of  experiments. 
Fr  =  0.43,  dry  transom: 

1.  -  Stern-quartering  view. 

2.  -  Stern  view. 

3.  -  Starboard  view. 

4.  -  Stern-quartering  view  compared  to  snapshot  of  experiments. 

5.  -  Stern  view  looking  down  compared  to  snapshot  of  experiments. 
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